Search for: All records

Abstract Volcanic evolution in ocean island settings is often controlled by variations in the chemistry and volumetric flux of magma from an underlying mantle plume. In locations such as Hawaiʻi or Réunion, this results in predictable variations in magma chemistry, the rate of volcanic activity, and the depth of magma storage with volcanic age and/or distance from the centre of plume upwelling. These systems, however, represent outliers in global plume volcanism due to their high buoyancy flux, frequent eruptions, and large distance from any plate boundary. Most mantle plumes display clear interaction with nearby plate boundaries, influencing the dynamics of solid plume material in the upper mantle and the distribution of melt across regions of active volcanism. Yet, the influence of plume–ridge interaction and plume–ridge distance on the structure, characteristics, and evolution of magma storage beneath ocean island volcanoes remains under constrained. In this study, we consider the evolution of magmatic systems in the Galápagos Archipelago, a region of mantle plume volcanism located 150–250 km south of the Galápagos Spreading Centre (GSC), focusing on the depth of magma storage during the eastward transport of volcanic systems away from the centre of plume upwelling. Geochemical analysis of gabbro xenoliths from Isla Floreana in the southeastern Galápagos suggest that they formed at ~2–2.5 Ma, when the island was located close to the centre of plume upwelling. These nodules, therefore, provide rare insights into the evolution of volcanic systems in the Galápagos Archipelago, tracking variations in the magma system architecture as the Nazca plate carried Isla Floreana eastwards, away from the plume centre. Mineral thermobarometry, thermodynamic modelling, and CO2 fluid inclusion barometry reveal that Isla Floreana’s plume-proximal stage of volcanic activity—recorded in the gabbro xenoliths—was characterized by the presence of high-pressure magma storage (>25 km), below the base of the crust. In fact, we find no petrological evidence that sustained, crustal-level magma storage ever occurred beneath Isla Floreana. Our results contrast with the characteristics of volcanic systems in the western Galápagos above the current centre of plume upwelling, where mid-crust magma storage has been identified. We propose that this change in magmatic architecture of plume-proximal volcanic centres in the Galápagos—from high-pressure mantle storage at 2.5 Ma to mid-crustal storage at the present day—is controlled by the variations in plume–ridge distance. Owing to the northward migration of the GSC, the distance separating the plume stem and GSC is not constant, and was likely <100 km at 2.5 Ma, significantly less than the current plume–ridge distance of 150–250 km. We propose that smaller plume–ridge distances result in greater diversion of plume-material to the GSC, ‘starving’ the eastern Galápagos islands of magma during their initial formation and restricting the ability for these systems to develop long-lived crustal magma reservoirs. 

Abstract Rapid-response petrological monitoring is a major advance for volcano observatories, allowing them to build and validate models of plumbing systems that supply eruptions in near-real time. The depth of magma storage has recently been identified as high-priority information for volcanic observatories, yet this information is not currently obtainable via petrological monitoring methods on timescales relevant to eruption response. Fluid inclusion barometry (using micro-thermometry or Raman spectroscopy) is a well-established petrological method to estimate magma storage depths and has been proposed to have potential as a rapid-response monitoring tool, although this has not been formally demonstrated. To address this deficiency, we performed a near-real-time rapid-response simulation for the September 2023 eruption of Kīlauea, Hawaiʻi. We show that Raman-based fluid inclusion barometry can robustly determine reservoir depths within a day of receiving samples—a transformative timescale that has not previously been achieved by petrological methods. Fluid inclusion barometry using micro-thermometric techniques has typically been limited to systems with relatively deep magma storage (>0.4 g/cm3 i.e.  > 7 km) where measurements of CO2 density are easy and accurate because the CO2 fluid homogenizes into the liquid phase. Improvements of the accuracy of Raman spectroscopy measurements of fluids with low CO2 density over the past couple of decades has enabled measurements of fluid inclusions from shallower magmatic systems. However, one caveat of examining shallower systems is that the fraction of H2O in the fluid may be too high to reliably convert CO2 density to pressure. To test the global applicability of rapid response fluid inclusion barometry, we compiled a global melt inclusion dataset (>4000 samples) and calculate the fluid composition at the point of vapor saturation ($${\mathrm{X}}_{{\mathrm{H}}_2\mathrm{O}}$$). We show that fluid inclusions in crystal hosts from mafic compositions (<57 wt. % SiO2)—likely representative of magmas recharging many volcanic systems worldwide—trap fluids with $${\mathrm{X}}_{{\mathrm{H}}_2\mathrm{O}}$$ low enough to make fluid inclusion barometry useful at many of the world’s most active and hazardous mafic volcanic systems (e.g. Iceland, Hawaiʻi, Galápagos Islands, East African Rift, Réunion, Canary Islands, Azores, Cabo Verde). 

We present DiadFit—an open-source Python3 tool for efficient processing of Raman spectroscopy data collected from fluid inclusions, melt inclusions and silicate melts. DiadFit is optimized to fit the characteristic peaks from CO2 fluids (Fermi diads, hot bands, 13C), gas species such as SO2, N2, solid precipitates (e.g. carbonates), and Ne emission lines with easily tweakable background positions and peak shapes. DiadFit's peak fitting functions are used as part of a number of workflows optimized for quantification of CO2 in melt inclusion vapour bubbles and fluid inclusions. DiadFit can also convert between temperature, pressure and density using various CO2 and CO2-H2O equations of state (EOS), allowing calculation of fluid inclusion pressures (and depths in the crust), conversion of homogenization temperatures from microthermometry to CO2 density, and propagation of uncertainties associated with EOS calculations using Monte Carlo methods. There are also functions to quantify the area ratio of the silicate vs. H2O region of spectra collected on silicate glasses to determine H2O contents in glasses and melt inclusions. Documentation and worked examples are available (https://bit.ly/DiadFitRTD, https://bit.ly/DiadFitYouTube). 

Abstract The iconic volcanoes of the Cascade arc stretch from Lassen Volcanic Center in northern California, through Oregon and Washington, to the Garibaldi Volcanic Belt in British Columbia. Recent studies have reviewed differences in the distribution and eruptive volumes of vents, as well as variations in geochemical compositions and heat flux along strike (amongst other characteristics). We investigate whether these along‐arc trends manifest as variations in magma storage conditions. We compile available constraints on magma storage depths from InSAR, geodetics, seismic inversions, and magnetotellurics for each major edifice and compare these to melt inclusion saturation pressures, pressures calculated using mineral‐only barometers, and constraints from experimental petrology. The availability of magma storage depth estimates varies greatly along the arc, with abundant geochemical and geophysical data available for some systems (e.g., Lassen Volcanic Center, Mount St. Helens) and very limited data available for other volcanoes, including many which are classified as “very high threat” by the USGS (e.g., Glacier Peak, Mount Baker, Mount Hood, Three Sisters). Acknowledging the limitations of data availability and the large uncertainties associated with certain methods, available data are indicative of magma storage within the upper 15 km of the crust (∼2 ± 2 kbar) beneath the main edifices. These findings are consistent with previous work recognizing barometric estimates cluster within the upper crust in many arcs worldwide. There are no clear offsets in magma storage between arc segments that are in extension, transtension or compression, although substantially more petrological work is needed for fine scale evaluation of storage pressures. 

The chemistry of erupted clinopyroxene crystals (±equilibrium liquids) have been widely used to deduce the pressures and temperatures of magma storage in volcanic arcs. However, the large number of different equations parameterizing the relationship between mineral and melt compositions and intensive variables such as pressure and temperature yield vastly different results, with implications for our interpretation of magma storage conditions. We use a new test dataset composed of the average Clinopyroxene-Liquid (Cpx-Liq) compositions from N = 543 variably hydrous experiments at crustal conditions (1 bar to 17 kbar) to assess the performance of different thermobarometers and identify the most accurate and precise expressions for application to subduction zone magmas. First, we assess different equilibrium tests, finding that comparing the measured and predicted Enstatite-Ferrosillite and KD (using Fet in both phases) are the most useful tests in arc magmas, whereas CaTs, CaTi and Jd tests have limited utility. We then apply further quality filters based on cation sums (3.95–4.05), number of analyses (N > 5) and the presence of reported H2O data in the quenched experimental glass (hereafter ‘liquid’) to obtain a filtered dataset (N = 214). We use this filtered dataset to compare calculated versus experimental pressures and temperatures for different combinations of thermobarometers. A number of Cpx-Liq thermometers perform very well when liquid H2O contents are known, although the Cpx composition contributes little to the calculated temperature relative to the liquid composition. Most Cpx-only thermometers perform very badly, greatly overestimating temperatures for hydrous experiments. These two findings demonstrate that the Cpx chemistry alone holds very little temperature information in hydrous systems. Most Cpx-Liq and Cpx-only barometers show similar performance to one another (mostly yielding root mean square errors [RMSEs] of 2–3.5 kbar), although the best Cpx-only barometers currently outperform the best Cpx-Liq barometers. We also assess the sensitivity of different equations to melt H2O contents, which are poorly constrained in many natural systems. Overall, this work demonstrates that Cpx-based barometry on individual Cpx only provides sufficient resolution to distinguish broad storage regions in continental arcs (e.g. upper, mid, lower crust). Significant averaging of Cpx compositions from experiments reported at similar pressures can reduce RMSEs to ~1.3–1.9 kbar. We hope our findings motivate the substantial amount of experimental and analytical work that is required to obtain precise and accurate estimates of magma storage depths from Cpx ± Liq equilibrium in volcanic arcs. 

Raman spectroscopy has become the tool of choice for analyzing fluid inclusions and melt inclusion (MI) vapor bubbles as it allows the density of CO2-rich fluids to be quantified. Measurements are often made at ambient temperature (Tamb ~18-25 °C), resulting in reported bulk densities between 0.2 and 0.7 g/mL despite that single-phase CO2 under these conditions is thermodynamically unstable and instead consists of a liquid (~0.7 g/mL), and a vapor phase (~0.2 g/mL). Here, we present results from experiments conducted at Tamb and 37 °C (above the CO2 critical temperature) on 14 natural CO2-rich MI bubbles from Mount Morning, Antarctica. Here, we show that at Tamb, laser power strongly affects the CO2 Raman spectrum of MI bubbles with bulk densities within the miscibility gap. High-power laser heating and low spectral resolution explain why published measurements have reported such bulk densities at Tamb even when using an instrument-specific calibration. 

Abstract The composition of clinopyroxene and clinopyroxene-liquid (Cpx-Liq) pairs are frequently used to calculate crystallization/equilibration pressures in igneous systems. While canonical uncertainties are often assigned to calculated pressures based on fits to calibration or test datasets, the sources of these uncertainties (and thus ways to reduce them) have not been rigorously assessed. We show that considerable uncertainties in calculated pressures arise from analytical error associated with Electron Probe Microanalyser (EPMA) measurements of Cpx. Specifically, low X-ray counts during analysis of elements with concentrations <1 wt% resulting from insufficient count times and/or low beam currents yield highly imprecise measurements (1σ errors of 10–40% for Na2O). Low analytical precision propagates into the calculation of pressure-sensitive mineral components such as jadeite. Using Monte Carlo approaches, we demonstrate that elemental variation resulting from analytical precision alone generates pressures spanning ~4 kbar (~15 km) for a single Cpx and ~6 kbar for a single Cpx-Liq pair using popular barometry expressions. In addition, analytical uncertainties in mineral compositions produce highly correlated arrays between pressure and temperature that have been previously attributed to transcrustal magma storage. Before invoking such geological interpretations, a more mundane origin from analytical imprecision must be ruled out. Most importantly, low analytical precision does not just affect the application of barometers to natural systems; it has also affected characterization of Cpx in experimental products used to calibrate and test barometers. The impact of poor precision on each individual measurement is often magnified by the small number of measurements made within experimental charges, meaning that low analytical precision and true variability in mineral compositions have not been sufficiently mediated by averaging multiple EPMA analyses. We compile the number of Cpx measurements performed in N = 307 experiments used to calibrate existing barometers, and N = 490 new experiments, finding ~45% of experiment charges were characterized by ≤5 individual Cpx analyses. Insufficient characterization of the true composition of experimental phases likely accounts for the fact that all Cpx-based barometers exhibit large errors (± 3 kbar) when tested using global experimental datasets. We suggest specific changes to analytical and experimental protocols, such as increased count times and/or higher beam currents when measuring low concentration elements in relatively beam resistant Cpx in experiments and natural samples. We also advocate for increasing the number of analyses per experimental charge, resolving interlaboratory analytical offsets and improving data reporting. Implementing these changes is essential to produce a more robust dataset to calibrate and test the next generation of more precise and accurate Cpx-based barometers. In turn, this will enable more rigorous investigation of magma storage geometries in a variety of tectonic settings (e.g. distinguishing true transcrustal storage vs. storage in discrete reservoirs). 

Abstract Magmas with matrix glass compositions ranging from basalt to dacite erupted from a series of 24 fissures in the first 2 weeks of the 2018 Lower East Rift Zone (LERZ) eruption of Kīlauea Volcano. Eruption styles ranged from low spattering and fountaining to strombolian activity. Major element trajectories in matrix glasses and melt inclusions hosted by olivine, pyroxene and plagioclase are consistent with variable amounts of fractional crystallization, with incompatible elements (e.g., Cl, F, and H₂O) becoming enriched by 4–5 times as melt MgO contents evolve from 6 to 0.5 wt%. The high viscosity and high H₂O contents (∼2 wt%) of the dacitic melts erupting at Fissure 17 account for the explosive Strombolian behavior exhibited by this fissure, in contrast to the low fountaining and spattering observed at fissures erupting basaltic to basaltic‐andesite melts. Saturation pressures calculated from melt inclusion CO₂‐H₂O contents indicate that the magma reservoir(s) supplying these fissures was located at ∼2–3 km depth, which is in agreement with the depth of a dacitic magma body intercepted during drilling in 2005 (∼2.5 km) and a seismically imaged lowVp/Vsanomaly (∼2 km depth). Nb/Y ratios in erupted products are similar to lavas erupted between 1955 and 1960, indicating that melts were stored and underwent variable amounts of crystallization in the LERZ for >60 years before being remobilized by a dike intrusion in 2018. We demonstrate that extensive fractional crystallization generates viscous and volatile‐rich magma with potential for hazardous explosive eruptions, which may be lurking undetected at many ocean island volcanoes.